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ABSTRACT 

Hot Jupiters are expected to be dark from both observations (albedo upper limits) and theory (alkali 
metals and/or TiO and VO absorption). However, only a handful of hot Jupiters have been observed 
with high enough photometric precision at visible wavelengths to investigate these expectations. The 
NASA Kepler mission provides a means to widen the sample and to assess the extent to which hot 
Jupiter albedos are low. We present a global analysis of Kepler-7b based on Q0-Q4 data, published 
radial velocities, and asteroseismology constraints. We measure an occultation depth in the Kepler 
bandpass of 44±5 ppm. If directly related to the albedo, this translates to a Kepler geometric albedo 
of 0.32±0.03, the most precise value measured so far for an exoplanet. We also characterize the 
planetary orbital phase lightcurve with an amplitude of 42±4 ppm. Using atmospheric models, we 
find it unlikely that the high albedo is due to a dominant thermal component and propose two solutions 
to explain the observed planetary flux. Firstly, we interpret the Kepler-7b albedo as resulting from 
an excess reflection over what can be explained solely by Rayleigh scattering, along with a nominal 
thermal component. This excess reflection might indicate the presence of a cloud or haze layer in 
the atmosphere, motivating new modeling and observational efforts. Alternatively, the albedo can be 
explained by Rayleigh scattering alone if Na and K are depleted in the atmosphere by a factor of 
10-100 below solar abundances. 

Subject headings: planetary systems - stars: individual (Kcpler-7, KIC 5780885, 2MASS 
19141956+4105233) - techniques: photometric 



1. INTRODUCTION 

More than 30 hot Jupiters benefit from observations 
of their emitted radiation from near to mid-infrared, 
where the measurement of their thermal emission is 
the most favorable. Spitzer made a significant con- 
tribution by producing a flurry of results allowing us 
to derive general proper ties of hot-Jupiter atmospheres 
(|Deming k. Sea ger 2009). Those planets are strongly 
irradiated by their host stars and their equilibrium 
temperatures were early estimated to be above 1000K 
l|Seager &; Sasselovll2000f h Those observations confirm 
that hot Jupiters efficiently reprocess the incident stellar 
flux into thermal reemission, exhibit i ng low flux at vis- 
ible wavelengths (IMarlev et al.|[l999HSeager et alj[2000t 
iSudarskv et al.ll2003t ). 

Characterization of transiting hot- Jupiter reflected 
light suffers from the scarcity of observations. The plan- 
etary to stellar flux ratio is of the order of 10~ 5 in the 
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visible, two orders of magnitude less than mid-infrared 
signatures. 

To date, eleven planets have an u pper limit constraint 
on th e ir geometric albedo : rBoo b (iCharbonneau et al.l 
Il999t iLeigh et alj I2003at iRodler et all |201C|L uAndb 



(Collier Came ron et alj|200l3) , HD75289Ab (Lei gh et al 
I2003bt IRodler et all I200lk HD2 09458b faowe et al 



Alonso et al 



Alonso et al 



2009; 



2010: 



2008) , CoRoT-lb 

2009) , CoRoT-2b 

B . HAT-P-7b (iChristiansen et alJ 12010 
20iqi . Kepler-5b ^Kipping fe Bakos ] 2011 1; 
2011 ). Kepler-6b jKipping fc Bakosl I20li1 : 
201l|), Ke pler-7b (IKipping & Bakosl 



Snellen et al 



Snell en et al 
iWelsh et a! 

iDesert et al 



2U1U ). Ke pler- yd ( Kipping sua 
HD189733b (|Berdvugina et al.N2011h . 



Desert et al 
20ll and 



Eight of them corroborate early theoretical predic- 
tions: with A g < 0.3 (3a upp er limit) hot Jupiters 
are da rk in the visible. However iCollier Cameron et al.l 
(2002) determined A g < 0.42 (3cr) f rom spectroscopy 
in the 380-650nm range for vAnd b, Kipping & Bakos 
(IMl reported A a = 0.38 ± 0.12^ for Kepler-7b and 
IBerdvugina et all (|2011[ ) determined a V-band albedo of 
A g = 0.28 ± 0.16 for HD189733b from polarimetry, sug- 
gesting dominance of reflected light over thermal emis- 
sion. 

Solar system giant planets have geometric albedos of 
0.32 (Ur anus) to 0.50 (Jup iter) in a bandpass similar to 
Kepler's (jKarko schka 1994). Those objects harbor bright 
cloud decks made of ammonia and water ice that are 
highly reflective at visible wavelengths. In contrast to the 
solar system giant planets, atmosphere models show that 

9 Determined in the Ke pler bandpass, wh ich has a >5% response 
between 423 and 897 nm IIKoch et ai.l feOlO') 
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the presence of alkali metals in hot- Jupiter atmospheres 
(Na and K) as well as TiO and VO (at the hotter range) 
causes significant absorption at visible wavelengths. 

We report in this letter the character ization of the hot 
Jupiter Kepler-7b (Lat ham et al.l[20Tol ). based on Kepler 
Q0-Q4 data. We present the photometry and data anal- 
ysis in Sect. 2, while corresponding results are shown in 
Sect. 3. Discussion and atmospheric analysis are finally 
presented in Sect. 4. 

2. OBSERVATIONS AND DATA ANALYSIS 

2.1. Kepler photometry 

Kepler-7b belongs to the first set of new planets pub- 
lished by the Kepler science team in early 2010. Kepler- 
7 b is a 4.8-day period hot Jupiter orbiting a T^=13.04 
sub-giant G sta r with M+ = 1 347 ± O.O7M and R* = 
1.843 ±0.07i? Q (|Latham et alJl2010j ). Like all objects lo- 
cated in the Kepler field, Kepler-7 benefits from nearly 
continuous photometric monitoring since mid-2009. 

We base our analysis on the Q0-Q4 quarters, which 
represent nearly one year of observations. Data recorded 
during each quarter are differentiated in short- and long- 
cadence timeseries, that are binnings per 58.84876s and 
29.4244min re spectively of the sam e CCD readouts. Five 
long-cadence (iJenkins et al.ll2010D and six short cadence 
(Gill iland et all 12010ft datasets are used as part of this 
study, representing 272,719 photometric datapoints and 
311.68 effective days of observations, out of which 175.37 
days have also been recorded in short cadence. We used 
the raw photometry for our purpose. 

Kepler-7 is a photometrically quiet star: apart from the 
4.88-day period transit signals, no evidence of significant 
stellar variability is apparent in the data. 

2.2. Data analysis 

For the purpose of this global analysis, we used 
the implementation of the Markov Chain Mont e- Carlo 
(MCMC) algorithm presented in iGillon et all (|2009L 
2010). MCMC is a Bayesian inference method based on 
stochastic simulations that samples the posterior prob- 
ability distributions of adjusted parameters for a given 
model. Our MCMC implementation uses the Metropolis- 
Hastings algorithm to perform this sampling. Our nom- 
inal model is based on a star and a transiting planet on 
a Keplerian orbit about their center of mass. 

Our global analysis was performed using 199 transit 
and occultation lightcurves in total, out of which 70 
were acquired in short cadence mode. We discarded 13 
lightcurves because of discontinuities due to spacecraft 
roll, change of focus, pointing offsets or safe mode events. 
Input data also include the 9 radial velocity points 
obtained from NOT/FIES ( Fiber-fed Eche l le Spe ctro- 
graph) that were published in lLatham et al.l (|2010l ). 

As the focus of this study is on using the transits 
and occultations to refine the system parameters, for the 
model fitting we use only the photometry near the eclipse 
events. Windows of width 0.6 days (12.3% of the or- 
bit) surrounding eclipses were used to measure the local 
out-of-transit baseline, while minimizing the computa- 
tion time. A dilution of 2.7±0.5% was determined from 
MMT/ARIES3 (ARizona Infrared imager and Echelle 

10 MMT is a joint facility of the Smithsonian Institution and 
University of Arizona. 



Spectrograph) observations and applied to both transit 
and occultation photometry. 

The excellent sampling of the transit lightcurve moti- 
vated us to fit for the limb-darkening (LD) coefficients. 
For this purpose, we assumed a quadratic law and used 
ci = 2u± + U2 and ci — u\ — 2u2 as jump parameters, 
where u\ and U2 are the quadratic coefficients. 

The MCMC has the following set of jump parameters: 
the planet/star flux ratio, the impact parameter b, the 
transit duration from first to fourth contact, the time of 
minimum light, the orbital period, K' = K\J\ — e 2 P 1 / 3 , 
where K is the radial-velocity semi-amplitude, the occul- 
tation depth, the two LD combinations c\ and C2 and the 
two parameters i/ecosw and -^/esinu;. A uniform prior 
distribution is assumed for all jump parameters. 

2.2.1. Model and systematics 

The tr ansit and occult a tion p hotometry are modeled 
with the iMandel fc Agoll (|2002f) model, multiplied by a 
second order polynomial accounting for stellar and in- 
strumental variability. Baseline model coefficients are 
determined for ea ch lightcurve with the SVD method 
(|Press et al.lll992H at each step of the MCMC. Co r relate d 
noise was a c count ed for following I Winn et"aT1 (|2008f k 
IGillon et al.l ((2010), to ensure reliable error bars on the 
fitted parameters. For this purpose, we compute a scal- 
ing factor based on the standard deviation of the binned 
residuals for each lightcurve with different time bins. The 
error bars are then multiplied by this scaling factor. We 
obtained a mean scaling factor of 1.04 for all photometry, 
denoting a negligible contribution from correlated noise. 
The mean global photometric RMS per 30-min bin is 96 
parts per million (ppm). 

2.2.2. Asteroseismology 

The data series for Kepler-7 contains 9 months of 
data at a cadence of 1 minute. The power spectrum 
shows a clear excess of power near 1.05mHz. The 
asteroseismic analysis of the data was performed us- 
ing the pipeline developed at the Kepler Asteroseis- 
mic Science Operations Cent e r as describ e d in detail by 
iChristensen-Dalsgaard et al.l (120081 I2010D : iHuber et all 
(|2009ft : iGilliland et al.l (j2^Tlh . Using the matched fil- 
ter approach we determine a value for the large separa- 
tion of 56/xHz. Locating the asymptotic frequency struc- 
ture in the folded power allows a robust identification of 
13 individual p-mode frequencies and an estimate of the 
scatter on those frequencies (0.9/iHz). The frequencies 
resulting from this analysis were fitted to stellar models 
in the same manner as in Christcnscn-Dalsgaar d~et al.l 
(2010), using also the effecti ve temperature and metal- 
licity ([Fe/H]) determined bv lLatham et al.1 (|2010l ). The 
models did not include diffusion and settling. Models 
were computed without overshoot from the convective 
core, as well as with overshoot of O.lflp and 0.2H p , where 
Hp is the pressure scale height at the edge of the convec- 
tive core. The observed frequencies and effective temper- 
ature were fitted to the models in a least-squares sense, 
resulting in a weighted average of the stellar properties. 
Interestingly, only models with overshoot provided ac- 
ceptable fits to the frequencies within the observed range 
of the effective temperature. 

We used the resulting stellar density (see Table 1) as 
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a Bayesian prior in the MCMC and the corresponding 
stellar mass to derive the system's physical parameters. 

2.2.3. Phase curve 

About 312 days of Kepler-7 observations are covered 
in Q0-Q4 data. This motivated us to search for the 
planetary phase signature. We first removed the tran- 
sits and fitted the long and short cadence data to re- 
move temporal linear trends. Stellar and instrumental in- 
duced modulation on the photometry was then removed 
by pre-whitening the raw data using Period04 software 
(|Lenz fe B rcger 2005). This step allowed us to filter out 
all frequencies below the orbital frequency and those that 
are not connected with the planetary orbit period. 

3. RESULTS 

We present Kepler-7 b's system parameters in Table 
1. Each value is the median of the marginal posterior 
distribution obtained for the relevant parameter. Error 
bars are the corresponding 68.3% probability interval. 
Figure 1 shows the phase-folded transit photometry of 
Kepler-7 b. 

Our results confirm the very low density of Kepler- 
7 b. The external constraint from asteroseismology mod- 
eling causes the planet radius to si gnificantly increase a s 
compared to the discovery paper (La tham et alj 12010). 
Our MCMC analysis yields a planetary radius of Rp = 
1.61 ± 0.02R Jup and a mass of M P = 0.44 ± 0.04M Jllp , 
which gives a surprisingly low density of pp = 0.14 ±0.01 
g cm -3 . Interes tingly, Kepler-7b's pr operties are close 
to WASP-17b's (jAnderson et all 120101) . the least dense 
transiting planet discovered so far with pp = 0.12±0.06 g 
cm -3 . Both planets are of similar mass and orbit evolved 
stars. 

We find a marginal orbital eccentricity of e = 0.001 ± 
0.001. 

We determine an occultation depth of 44±5 ppm. The 
corresponding phase-folded lightcurve is shown on bot- 
tom panel of Fig. 2. 

Finally we find an orbital phase curve of 42±4 ppm am- 
plitude that is consistent with the occultation depth and 
phased with Kepler-7 b's transits and occultations. We 
model the phase cur ve assuming a Lambert law phase- 
dependent flux-ratio ()Sobolevlll975l ) : 



(a) = A g 



R,. 



sin(a) + (ir — a) cos(a) 



where a is the orbital phase, A g is the geometric albedo, 
R p is the planetary radius and a is the orbital semi- 
major axis. A value of A g = 0.31 ± 0.03 is deduced 
from the phase curve using R p and a values from Ta- 
ble 1. Although a perfectly reflecting Lambertian sphere 
has A g = |, the Lambert law can represent a scatterer 
which takes on a lower A a val ue for an atmosphere with 
some absorption (jSeagcr 20lQ). The resulting lightcurve 
is shown on top panel of Fig. 2, with the best-fit model 
superimposed. No ellipsoidal variations are detected 
(— 1 ± 3ppm). 

4. THE HIGH ALBEDO OF KEPLER-7 B 

If directly related to the albedo, the measured occul- 
tation depth of 44±5 ppm translates to a 0.32±0.03 geo- 
metric albedo as measured by Kepler. The Kepler band- 
pass encompasses a large range of wavelengths, from 0.4 



to 0.9 /xm. Albedo values reported in this paper for 
Kepler-7 b are averaged over this spectral domain. For 
the highly irradiated exoplanets, a significant part of the 
thermal emission leaks into the red end of this bandpass, 
making the occultation depth larger. With a 4.9-d long 
orbit, Kepler-7 b would not be expected to be one of the 
hottest giant planets found to date. However, its host 
star with a T cff =5933K, 1.4 M and 2.0 R Q , is about 
4.5 times more luminous than the Sun. These compen- 
sating factors make it necessary to estimate the relative 
contributions of thermal emission and reflected light to 
the occultation depth. 

The possible relative contributions of thermal emis- 
sion and reflected light to the observed flux in the Ke- 
pler bandpass are shown in Fig. 3. The thermal emis- 
sion is represented by an effective brightness tempera- 
ture (Tb), and the reflected light by the geometric albedo 
(A g ). Also shown are the dayside equilibrium tempera- 
tures corresponding to atmosphere with efficient versus 
inefficient energy redistribution (dotted lines). The de- 
generacy between Tp and A g is evident. The observed 
occultation depth allows for geometric albedos as high 
as 0.35 for Tp = 1500K in the Kepler bandpass. On 
the other hand, allowing for a zero geometric albedo re- 
quires an extremely high Tp of 2500-2600K, which is ~ 
400K higher than the maximum equilibrium tempera- 
ture. Completely breaking the degeneracy between A g 
and Tp requires additional observations in the visible 
and near-infrared. However, tentative constraints can be 
placed on the various sources of opacity and scattering 
using a physical model atmosphere. 

The thermal and reflection spectra of the planetary 
atmosphere depend on the various atomic and molec- 
ular opacities, sources of scattering, and the temper- 
ature structure, all of which constitute a large num- 
ber of free parameters under-constrained by the single 
data point available. Nevertheless, the precise mea- 
surement allows us to constrain regions of the param- 
eter space, given plausible assumptions about energy 
balance and atmospheric chemistry. We use the exo- 
planet atmospheric mod e ling a nd retrieval technique of 
Madhusu dhan k, Seagerl (|2009). The model computes 
line-by-line radiative transfer in a plane-parallel atmo- 
sphere, with the assumption of hydrostatic equilibrium 
and global energy balance, and parametrized chemical 
composition and temperature structure. A key aspect 
of the model is the flexibility to explore a wide range of 
molecular abundances and pressure-temperature (P-T) 
profiles, without any assumptions of chemical or radia- 
tive equilibrium. We consider all the major sources of 
opacity in the visible and near-infrared: H2-H2 collision- 
induced absorption, Na, K, Rayleigh scattering, TiO and 
VO ( with a condensation curve), H 2 Q, CO, CH 4 and 
C0 2 ([Christiansen et all l20l0t IMadhusudhan fc Seagerl 
I2010I I2oTTIl ~ However, given the availability of the Ke- 
pler observation alone in the present work, the number 
of free parameters far exceed the single data point. Con- 
sequently, we explore a range of 1-D P-T profiles and 
chemical composition to explore possible explanations of 
the observed flux in the Kepler bandpass. 

Our results indicate that a high geometric albedo is 
the most plausible explanation of the observed eclipse 
depth. The major sources of opacity in the Kepler band- 
pass are atomic Na and K, molecular TiO and VO where 
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Fig. 2. — Top : Kepler-7b orbital phase curve with best-fit model (see text) superimposed. Transits are omitted. Bottom : Kepler-7b 
phase-folded occultation lightcurve with best-fit model. Binned per 15min. 
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TABLE 1 
Kepler- 7 system parameters 



Parameters 



Value 



Jump parameters 



Planet/star area ratio (R p /R s ) 2 

V = acosi/R^ [R+] 

Transit width [d] 

T - 2450000 [HJD] 

Orbital period P [d] 

RV K' [ms-M 1 / 3 ] 

Y^cosw 

v^sinw 

ci = 2«i + u 2 

C 2 = Ml - 2« 2 

Occultation depth 



006772 +0 '°°°° 18 
U.UUOf (4_ 000021 

U.OODO_ g 63 

n 01 777+0.00023 
U.Z1 ( ( / -0.00021 

4967.275991°;°°°^ 
4.8854830l°;°~ 

7 o 1+6.7 

n 0376+ a ° 143 
U - Ud ' CI -0.0153 

„+0.0026 
^-0.0022 
,+0.011 
-0.011 

-0.1431°;°^ 
0.0000441° = 



-0.0016] 
0.922 4 



Deduced stellar parameters 



■hi 
U2 



Density p 3ta r [g cm~ 3 ] 
Surface gravity log<?„ [cgs] 
Mass M* [M©] 
Radius [Rq] 



0.344] 
0.232^ 
0.23lj 
3.960^ 



1.36 
2.02 



+0.03 
-0.03 
+0.02 
0.02 



007 
007 
009 
009 
004 
004 
00.1 
005 



Asteroseismic parameters 



Density p sta r [g cm" 
Mass M star [Mq] 
Radius R 3ta r [Re] 
Stellar age [Gyr] 
Surface gravity logg. 



legsj 



0.252 ± 0.003 
1.359 ± 0.031 
1.966 ± 0.013 
3.3 ± 0.4 
3.984 ± 0.006 



Deduced planet parameters 



RV K [ms- 1 ] 

btransit [R*] 
^occultation [R*] 

T occuHation - 2450000 [HJD] 
Orbital semi-major axis a [AU] 
Orbital inclination i [deg] 
Orbital eccentricity e 
Argument of periastron ui [deg] 
Density pp [g enr 3 ] 
Surface gravity log gp [cgs] 
Mass M P [M Jup ] 
Radius Rp [Rj up ] 



43.1 



+3.9 



-4.0 
7 +0.006 
' -0.006 
-+0.006 
-0.006 

4974.6086l°;°°|° 
0.062461°°°°^ 

3+0.076 
-0.074 



0.557] 
0.556^ 



85.18] 



0.00ll°;°°i 
357.ll|:| 



0.14 



+0.01 



-0.01 

,+0.04 

2 -0.04 

3+0.041 

-0.042 

1 fi14+0- 015 
l.D14_o_ 018 



2.62] 
0.443] 



G 



0.4 r 




1600 1800 2000 2200 2400 2600 

Brightness Temperature T B [K] 

Fig. 3. — Model-independent domain of allowed geometric albedo 
A g (reflected light) and brightness temperature Tb (thermal emis- 
sion) combinations, constrained by the observed Kepler occultation 
depth. The shaded area represents the la confidence domain found 
using a Kurucz spectrum for the star and blackbody emission for 
the planet. The dotted lines depict the equilibrium temperature 
assuming full redistribution (left) and no redistribution (right). 



temperatures are higher than the condensation tempera- 
ture, Rayleigh scattering, and any possible contribution 
due to clouds and/or hazes. Consequently, we examine 
the constraints on each of these opacity sources due to 
the Kepler point. We assume all other molecules (H2O, 
CO, CH4 and CO2) to be in chemical equilibrium, assum- 
ing solar abundances. We find that the observed Kepler 
flux can be explained under three scenarios, shown in the 
three panels of Fig. 4. 

In a first scenario, we consider a model where all the 
major absorbers, i.e. Na, K, TiO, and VO, are in chem- 
ical equilibrium with solar abundances, and we nomi- 
nally vary only the P-T profile to find a close match to 
the data. In this scenario, we find that the observed 
flux cannot be accounted for by thermal emission and 
Rayleigh scattering alone, both of which together provide 
a net flux contrast of ~ 20 ppm, with Rayleigh scatter- 
ing contributing a geometric albedo of 0.15. In principle, 
one might expect that hotter temperature profiles might 
lead to greater thermal emission which could contribute 
to the observed flux. However, warmer P-T profiles in- 
tercept longer absorption columns of TiO /VO which fur- 
ther lower the emergent thermal flux. Therefore, under 
the assumption of solar abundances of all species, excess 
flux in the form of reflected light, potentially from clouds 
and/or hazes, is required to explain the observed Kepler 
flux, implying a net geometric albedo of ~ 0.3. 

In a second scenario, we investigate if thermal emission 
can be a predominant contributor to the observed Kepler 
flux. The strongest absorber in the redder regions of the 
Kepler bandpass, where thermal emission dominates, is 
TiO, followed by VO. As discussed in the equilibrium 
scenario above, TiO absorption precludes high bright- 
ness temperatures due to thermal emission in the Kepler 
bandpass. Conversely, we find that thermal emission can 
contribute substantially to the observed flux, if TiO and 
VO are assumed to be depleted by over 10 3 in the lower 
atmosphere (P ~ 0.1 — 1 bar). Such a scenario just man- 
ages to fit the data at the lower error bar, with equal 



contributions from thermal emission and Rayleigh scat- 
tering (with a geometric albedo of 0.17). However, two 
problems confront this scenario. Firstly, requiring such 
a high thermal flux implies that all the stellar incident 
flux on the dayside of the planet must be reradiated on 
the same side, with almost no energy recirculation to the 
night side. The Kepler phase curve of Kepler-7b, shown 
in Fig. 2, could support this scenario if all the flux were 
purely due to thermal emission. However, because there 
is still a >50% contribution due to Rayleigh scattering 
required to explain the net flux at occultation, the Ke- 
pler phase curve cannot be interpreted as solely due to 
a large day-night temperature contrast. Furthermore, 
at the high pressures (P ~ 1 bar) probed by the Ke- 
pler bandpass, the temperature should be homogenized 
across day and night side, exhibiting small thermal or- 
bital phase variation. Secondly, at the high tempera- 
ture of ~2600K probed by Kepler, it seems unlikely that 
TiO/VO can be depleted by factors of 10 3 below chemical 
equilibrium at ~ 1 bar pressure, although the difficulty 
of sustaining TiO/VO at high altitudes (P <-10~ 3 bar) 
has been reported in literature (Spiegel ct al. 2009). 

In our final scenario, we investigate if Rayleigh scat- 
tering alone can contribute dominantly to the observed 
flux. The strongest absorbers of Rayleigh scattered light, 
as shown in the first two panels of Fig. 4, are atomic Na 
and K. We find that if we allow depletion of Na and K 
by a factor of 10-100 of the equilibrium composition, the 
observed flux in the Kepler bandpass can be fit extremely 
well with a geometric albedo of 0.32 from Rayleigh scat- 
tering alone along with a nominal contribution due to 
thermal emission. The resultant model also yields effi- 
cient day-night energy circulation. 

We finally note that the large planetary radius and 
very low density would make appealing the hypothe- 
sis of Kepler-7b being a very young planet that would 
still be in its cooling phase. However, the asteroseismol- 
ogy results presented in Sect. 2 produce a stellar age of 
3.3±0.4Gyr, which argues for a pl anetary evolutionar y 
state beyond the collapsing phase (jFortnev et al J 12005( 1. 
and a negligible contribution from internal heat to the 
occultation depth. 

5. CONCLUSIONS 

Given the three scenarios described in Sect. 4, we in- 
terpret the Kepler observed planetary flux as due to a 
combination of Rayleigh sc attering and the presence of 
clouds or a haze layer (e.g., iLecavelier Pes Etangs et all 
120081 : iSing et all 120091 ) in the atmosphere of Kepler-7b, 
yielding an averaged geometric albedo of ~ 0.3 in the 
Kepler bandpass. A detailed cloud or haze model is be- 
yond the scope of the present work. Our results motivate 
new modeling and observational efforts to investigate the 
nature of clouds and hazes that might be possible in a 
low gravity atmosphere such as that of Kepler-7 b. 

Authors thank the Kepler Giant Planet Working 
Group, P.-O. Quirion, M. Holman, D. Ragozzine, J. 
Jenkins, J.-M. Desert, B. Benneke and D. Latham for 
useful discussions. We thank E. Dunham, W. Cochran 
and the referee for helpful comments that improved the 
manuscript. B.-O.D. acknowledges support from the 
Swiss NSF and thanks kindly R. Stewart, P.L. Vidale 



7 




Fig. 4. — Model spectra of Kepler-7b in the Kepler bandpass. Each panel corresponds to each scenario (1 to 3, left to right) described in 
Sect. 4. The Kepler data point is shown in red. In each panel, the black solid curve shows the net emergent spectrum. The net emergent 
flux integrated in the Kepler bandpass is shown in green. The cyan and purple curves show the contributions of thermal emission and 
reflected light to the emergent spectrum. The black dashed line shows a black body spectrum of the planet at 2550K, divided by a Kurucz 
model spectrum for the star. The relevant Pressure- Temperature profiles are shown for each atmospheric model in the insets, with the TiO 
condensation curve in dotted line. 
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